Searching for a Stochastic Background of Gravitational Waves 

with LIGO 



B. Abbott^^ R. Abbott^^ R. Adhikari^^ J. Agresti^^ R Ajith^, B. Allen^^ R. Amini^ 
S. B. Anderson^^, W. G. Anderson*^^, M. Araya^^, H. Armandula^^, M. Ashley^, S Aston^'^, 
C. AulbertS S. Babak\ S. Ballmer^^, B. C. Barish^^, C. Barker^^ D. Barker^^ B. Barr^e, 
^ i P. Barriga"^^, M. A. Barton^^, K. Bayer^^, K. Belczynski^^ , J. Betzwieser^^, P. Beyersdorf^^, 

O : B. Bhawal^^^ I. A. Bilenko^^ G. Billingsley^^^ E. Black^^ K. Blackburn^^^ L. Blackburn^^^ 

<^ : D. Blair4\ B. Blandl^ L. Boguel^ R. Bork^^ S. Bose^^ P. R. Brady^^^ V. B. Braginsky^^ 

^1 J. E. Brau^^, A. Brooks^^, D. A. Brown^^, A. Bullington^^, A. Bunkowski^, A. Buonanno^^ 
C/^ : R. Burman^\ D. Busby^^^ R. L. Byer^^, L. Cadonati^^^ G. CagnolP, J. B. Camp^o, 

J. Cannizzo^o, K. Cannon^^^ C. A. Cantley^^, J. Cao^^^ L. Cardenas^^^ M. M. Casey^^, 
C. Cepeda^^^ P. Charltoni^, S. Chatterji^^ S. Chelkowski^, Y. Chen\ D. Chinas, E. Chin^^ 
<^ : J. Chow^, N. Christensen^ T. Cokelaer^ C. N. Colacino^^, R. ColdwelF, D. Cookl^ 

\D ; T. Corbitt^^ D. Coward^\ D. Coyne^^ J. D. E. Creighton^^ T. D. Creighton^^ 

D. R. M. Crooks^s, A. M. Cruise^^, A. Cumming^e, C. Cutler^ J. Dalrymple^^, 



(N 



O 



00 I E. D'Ambrosio^^, K. Danzmann^^'^, G. Davies^, G. de Vine^, D. DeBra^^, J. Degallaix"^^ , 

§ ; V. Dergachev^s, S. DesaP, R. DeSalvo^^^ S. Dhurandar", A. Di Credico^^ M. Diazes, 

O ; J. Dickson^, G. Diederichs^S A. Dietz^^, E. E. Doomes^^ R. W. P. Drever^ J.-C. Dumas^^ 
-g^; R. J. Dupuis^^^ P. Ehrens^^^ E. Elliffe^^, T. Etzel^^^ M. Evans^^^ T. Evansl^ S. Fairhurst^^ 

I ; Y. Fan^\ M. M. Fejer^^, L. S. Finn^s, N. Fotopoulos^^ A. Franzen3\ K. Y. Franzen3^ 

^ ; R. E. Frey39, T. Fricke^o, P. Fritschel^^^ V. V. Frolovl^ M. FySe^^, J. Garofolil^ 

g ; 1. Gholami^ J. A. Giaime^^ S. Giampanis^°, K. Goda^^^ E. Goetz^^, L. Goggin^^ 

G. Gonzalezes, S. Gossler^, A. Grant^^, S. Gras^^, C. Grayl^ M. Gray^, J. Greenhalgh^^, 
^ ^ A. M. Gretarsson^, D. Grimmett^^, R. Grosso^^, H. Grote^, S. Grunewald^, M. Guenther^^, 
^ [ R. Gustafson^^, B. Hage^^, C. Hanna^^, J. Hanson^^, C. Hardham^^, J. Harms^, G. Harry^^ 
E. Harstad^^ T. Hayler^^, J. Heefner^^ 1. S. Heng^e, A. HeptonstalP, M. Heurs^^ 
M. Hewitson^, S. Hild^i, N. Hindmanl^ E. Hirose^^ D. Hoak^^ P. Hoangi^, D. Hosken^^^ 
J. Hough^e, E. Howell^S D. Hoyland^^, W. Hua^^, S. Huttner^^, D. Ingrami^ M. Ito^^, 
Y. Itoh^^ A. Ivanovi2, D. Jackrel^^, B. Johnson^^ W. W. Johnsoni^, D. 1. Jones^^, 
G. Jones^, R. Jones^^, L. Ju'^^ P. Kalmus^, V. Kalogera^^, D. Kasprzyk^^, 
E. Katsavounidis^^, K. Kawabe^^, S. Kawamura^^, F. Kawazoe^^, W. Kells^^, 
F. Ya. Khalilil^ A. Khanl^ C. Kim^^, P. King^^, S. Klimenko^^, K. Kokeyama^^, 
V. Kondrashov^^, S. Koranda^^, D. Kozak^^, B. Krishnan^, P. Kwee'^^, P. K. Lam^, 
M. Landry^"^, B. Lantz^^, A. Lazzarini^^, B. Lee^^, M. Lei^^, V. Leonhardt^^, I. Leonor^^, 
K. Libbrecht^^, P. Lindquist^^, N. A. Lockerbie^"^, M. Lormand^^, M. Lubinski^^, 
H. Liick3i'2, B. Machenschalk^ M. Maclnnis^^ M. Mageswaran^^ K. Mailand^^ 
M. Malec^\ V. Mandic^^ S. Marka^ J. Markowitz^^ E. Maros^^ I. Martin^^ J. N. Marx^^ 



- 2 - 



K. Mason^^ L. Matone^ N. Mavalvala^^ R. McCarthy^^ D. E. McClelland^ 
S. C. McGuire^s, M. McHugli^^ K. McKenzie^, J. W. C. McNabb^s, T. Meier^^, 

A. Melissinos^", G. Mendell^^ R. A. Mercer^^ S. Meshkov^^ E. Messaritaki'^^^ 
C. J. Messenger^e, D. Meyers^^^ E. Mikhailoyi^ S. Mitra^^ V. R Mitrofanov^^ 

G. Mitselmakher3^ R. Mittleman^^ O. Miyakawa^^ S. Mohanty^^, G. Moreno^^, 

K. Mossavi^, C. MowLowry^ A. Moylan^ D. Mudge3^ G. Mueller^^ H. Miiller-Ebhardt^ 
S. Mukherjee^^, J. Munch^^, P. Murray^^, E. Myers^^, J. Myers^^, G. Newton^^, 
K. Numata^o, B. 0'Reillyl^ R. O'Shaughnessy^^, D. J. Ottawayl^ H. Overmieri^^ 
B. J. Owen^^ Y. Pan^ M. A. Papa^'^^ V. Parameshwaraiah^^, M. Pedraza^^ S. Penn^°, 
M. Pitkin^e, M. V. Plissi^^, R. Prix\ V. Quetscllke^^ F. Raab^'', D. Rabeling^ 
H. Radkins^^, R. Rahkola3^ M. Rakhmanov^s, K. Rawlins^^ S. Ray-Majumder^^ V. Re^^, 
H. Rehbein^, S. Reid^^, D. H. Reitze^^, L. Ribichini^, R. Riescn^^ K. Rilcs^s, B. Rivera^^ 
D. I. Robcrtson^^, N. A. Robertson^^'^^, C. Robinson^, S. Roddy^^, A. Rodriguez^^, 
A. M. Rogan^^ J. Rollins^ J. D. Romano^ J. Romiel^ R. Routc^^, S. Rowan^^, 

A. Riidiger^, L. Ructi'\ P. Russell^^ K. Ryan^\ S. Sakata^^, M. Samidi^^ 
L. Sancho de la Jordana^^ V. Sandbergl^ V. Sannibale^^^ S.Saraf^^ P. Sarin^^^ 

B. S. Sathyaprakash^, S. Sato2^ P. R. Saulson^^, R. Savagel^ S. Schediwy^^ R. Schilling^, 
R. SchnabeP, R. Scllo^ield=^^ B. F. Scllutzl'^ P. Scllwinbergl^ S. M. Scott^ S. E. Seader^^ 

A. C. Searle^, B. Sears^^ F. Seifert^, D. Sellers^^ A. S. Sengupta^ P. Shawhan^^ 
B. Sheard^, D. H. Shoemakeri^, A. Sibleyl^ X. Siemens^^ D. Siggl^ A. M. Sintes^^'i, 

B. Slagmolen^ J. Slutsky^^ J. Smith^ M. R. Smith^^ P. Sneddon^^ K. Somiya^'\ 

C. Speake'■^^ O. Spjeldl^ K. A. Strain^^, D. M. Stromas, A. Stuver^^, T. Summerscales2^ 
K. Sun26, M. Sungi^ P. J. Sutton^^ D. B. Tanner^s, M. Tarallo^^ R. Taylor^^ R. Taylor^^, 
J. Tllackerl^ K. A. Tllorne2^ K. S. Tllorne^ A. Thiiring^i, K. V. Tokmakov^^ C. Torres2^ 

C. Torrie^^^ G. Traylorl^ M. Trias3^ W. lyier^^ D. UgolinP, C. Ungarelli^^, 

H. Valllbrucll3^ M. Vallisneri^, M. Varvellal^ S. Vass^^^ A. Vecchio=^^ J. Veitch^e, 

P. Veitcll3^ S. Vigeland^ A. Villari^, C. Vorvickl^ S. P. Vyachanin^^, S. J. Waldmani^, 

L. Wallace^^ H. Ward^^, R. Ward^l K. Watts^^ D. Webbcr^^ A. Weidncr^, 
A. Wcinstcin^^, R. Wciss^^, S. Wen^^, K. Wette^, J. T. Whe\an^^'\ D. M. Whitbeck^^, 
S. E. Whitcombi2, B. F. Whiting^^, C. Wilkinsoni^ P. A. Willemsl^ B. Willke3i'2, 

I. Wilmut23, W. Winkler^, C. C. Wipf ^ S. Wise^s, A. G. Wiseman^^ G. Woan^e, 

D. Woods^^ R. Wooley^^ J. Worden^^ W. Wu^^, I. Yakushini^ H. Yamamoto^^ Z. Yan^\ 
S. Yosllida2^ N. Yunes2^ M. Zanolinl^ L. Zllang^^ C. Zliao^\ N. Zotov^^ M. Zucker^^ 

H. zur Miihlen^^ J. Zweizig^^, 



The LIGO Scientific Collaboration, http://www.ligo.org 



-3- 



^Albert-Einstein-Institut, Max-Planck-Institut fiir Gravitationsphysik, D-14476 Golm, Germany 

^Albcrt-Einstcin-Institut, Max-Planck-Institut fiir Gravitationsphysik, D-30167 Hannover, Germany 

■'Australian National University, Canberra, 0200, Australia 

^California Institute of Technology, Pasadena, CA 91125, USA 

^Caltech-CaRT, Pasadena, CA 91125, USA 

f^Cardiff University, Cardiff, CF2 3YB, United Kingdom 

^Carleton College, Northfield, MN 55057, USA 

^Columbia University, New York, NY 10027, USA 

^Embry-Riddle Aeronautical University, Prescott, AZ 86301 USA 
i°Hobart and William Smith Colleges, Geneva, NY 14456, USA 
^^Inter-University Centre for Astronomy and Astrophysics, Pune - 411007, India 
i^LIGO - Cahfornia Institute of Technology, Pasadena, CA 91125, USA 
i^LIGO - Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
i^LIGO Hanford Observatory, Richland, WA 99352, USA 
i^LIGO Livingston Observatory, Livingston, LA 70754, USA 
-•^^Louisiana State University, Baton Rouge, LA 70803, USA 
^"^Louisiana Tech University, Ruston, LA 71272, USA 
^^Loyola University, New Orleans, LA 70118, USA 
^^Moscow State University, Moscow, 119992, Russia 
20NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 
^^National Astronomical Observatory of Japan, Tokyo 181-8588, Japan 
^^Northwestern University, Evanston, IL 60208, USA 

^^Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX United Kingdom 
^^Southeastern Louisiana University, Hammond, LA 70402, USA 
^^Southern University and A&M College, Baton Rouge, LA 70813, USA 
^^Stanford University, Stanford, CA 94305, USA 
^^Syracuse University, Syracuse, NY 13244, USA 

^^The Pennsylvania State University, University Park, PA 16802, USA 

^^The University of Texas at Brownsville and Texas Southmost College, Brownsville, TX 78520, USA 
30Trinity University, San Antonio, TX 78212, USA 



-4- 



ABSTRACT 



The Laser Interferometer Gravitational- wave Observatory (LIGO) has per- 
formed the fourth science run, S4, with significantly improved interferometer sen- 
sitivities with respect to previous runs. Using data acquired during this science 
run, we place a limit on the amplitude of a stochastic background of gravitational 
waves. For a frequency independent spectrum, the new limit is ^Iqw < 6.5 x 10~^. 
This is currently the most sensitive result in the frequency range 51-150 Hz, with 
a factor of 13 improvement over the previous LIGO result. We discuss comple- 
mentarity of the new result with other constraints on a stochastic background of 
gravitational waves, and we investigate implications of the new result for different 
models of this background. 

Subject headings: gravitational waves 



1. Introduction 

A stochastic background of gravitational waves (GWs) is expected to arise as a super- 
position of a large number of unresolved sources, from different directions in the sky, and 
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with different polarizations. It is usually described in terms of the GW spectrum: 

i!Gw(/) = ^%^. (1) 

Pc d} 

where dpcw is the energy density of gravitational radiation contained in the frequency range 
f to f + df (Allen & Romano 1999), pc is the critical energy density of the Universe, and 
/ is frequency (for an alternative and equivalent definition of ^cwif) see, for example, 
(Baskaran et al. 2006)). In this paper, we will focus on power-law GW spectra: 

Here, Qa is the amplitude corresponding to the spectral index a. In particular, Qq denotes 
the amplitude of the frequency-independent GW spectrum. 

Many possible sources of stochastic GW background have been proposed and several 
experiments have searched for it (see (Maggiore 2000; Allen 1996) for reviews). Some of the 
proposed theoretical models are cosmological in nature, such as the amplification of quantum 
vacuum fiuctuations during infiation (Grishchuk 1975), (Grishchuk 1997), (Starobinsky 1979), 
pre-big-bang models (Gasperini & Veneziano 1993), (Gasperini & Veneziano 2003), (Buonanno et al. 1997 
phase transitions (Kosowsky et al. 1992), (Apreda et al. 2002), and cosmic strings (Caldwell & Allen 1992^ 
(Damour & Vilenkin 2000), (Damour & Vilenkin 2005). Others are astrophysical in nature, 
such as rotating neutron stars (Regimbau & de Preitas Pacheco 2001), supernovae (Coward et al. 2002) 
or low-mass X-ray binaries (Cooray 2004). 

A number of experiments have been used to constrain the spectrum of GW background 
at different frequencies. Currently, the most stringent constraints arise from large-angle cor- 
relations in the cosmic microwave background (CMB) (Allen & Koranda 1994; Turner 1997), 
from the arrival times of millisecond pulsar signals (Jcnet et al. 2006), from Dopplcr tracking 
of the Cassini spacecraft (Armstrong et al. 2003), and from resonant bar GW detectors, such 
as Explorer and Nautilus (Astone et al. 1999). An indirect bound can be placed on the total 
energy carried by gravitational waves at the time of the Big-Bang Nucleosynthesis (BBN) 
using the BBN model and observations (Kolb & Turner 1990; Maggiore 2000; Allen 1996). 
Similarly, (Smith et al. 2006a) used the CMB and matter spectra to constrain the total 
energy density of gravitational waves at the time of photon decoupling. 

Ground-based interferometer networks can directly measure the GW strain spectrum in 
the frequency band 10 Hz - few kHz, by searching for correlated signal beneath uncorrelated 
detector noise. LIGO has built three power-recycled Michelson interferometers, with a Fabry- 
Perot cavity in each orthogonal arm. They are located at two sites, Hanford, WA, and 
Livingston Parish, LA. There are two collocated interferometers at the WA site: HI, with 
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4km long arms, and H2, with 2km arms. The LA site contains LI, a 4km interferometer, 
similar in design to HI. The detector configuration and performance during LIGO's first 
science run (SI) was described in (Abbott et al. 2004a). The data acquired during that run 
was used to place an upper limit of Qq < 44.4 on the amplitude of a frequency independent 
GW spectrum, in the frequency band 40-314 Hz (Abbott ct al. 2004b). This limit, as well as 
the rest of this paper, assumes the present value of the Hubble parameter Hq = 72 km/s/Mpc 
or, equivalently, /iioo = i^o/(100 km/s/Mpc) = 0.72 (Bennet et al. 2003). The most recent 
bound on the amplitude of the frequency independent GW spectrum from LIGO is based on 
the science run S3: Qq < 8.4 x 10~^ for a frequency-independent spectrum in the 69-156 Hz 
band (Abbott et al. 2005). 

In this paper, we report much improved limits on the stochastic GW background around 
100 Hz, using the data acquired during the LIGO science run S4, which took place between 
February 22, 2005 and March 23, 2005. The sensitivity of the interferometers during S4, 
shown in Figure 1, was significantly better as compared to S3 (by a factor 10 x at certain 
frequencies), which leads to an order of magnitude improvement in the upper limit on the 
amplitude of the stochastic GW background: flo < 6.5 x 10~^ for a frequency-independent 
spectrum over the 51-150 Hz band. 

This limit is beginning to probe some models of the stochastic GW background. As 
examples, wc investigate the implications of this limit for cosmic strings models and for 
prc-big-bang models of the stochastic gravitational radiation. In both cases, the new LIGO 
result excludes parts of the parameter space of these models. 

The organization of this paper is as follows. In Section 2 we review the analysis procedure 
and present the results in Section 3. In Section 4, we discuss some of the implications of our 
results for models of a stochastic GW background, as well as the complementarity between 
LIGO and other experimental constraints on a stochastic GW background. We conclude 
with future prospects in Section 5. 



The cross-correlation method for searching for a stochastic GW background with pairs 
of ground-based interferometers is described in (Allen & Romano 1999). We define the fol- 
lowing cross-correlation estimator: 



2. Analysis 



2.1. 



Cross-Correlation Method 
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H1 
H2 




Frequency (Hz) 

Fig. 1. — Typical strain amplitude spectra of LIGO interferometers during the science run 
S4 (solid curves top-to-bottom at 70 Hz: H2, LI, HI). The black dashed curve is the LIGO 
sensitivity goal. The gray dashed curve is the strain amplitude spectrum corresponding to the 
limit presented in this paper for the frequency-independent GW spectrum flo < 6.5 x 10~^. 
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where 5t is a finite-time approximation to the Dirac delta function, §1 and §2 are the Fourier 
transforms of the strain time-series of two interferometers, and Q is a filter function. As- 
suming that the detector noise is Gaussian, stationary, uncorrelated between the two inter- 
ferometers, and much larger than the GW signal, the variance of the estimator Y is given 
by: 
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where Pi{f) are the one-sided power spectral densities (PSDs) of the two interferometers and 
T is the measurement time. Optimization of the signal-to-noise ratio leads to the following 
form of the optimal filter (Allen & Romano 1999): 

n(f\ Kf 7(/)'S'gw(/) /p-v 
««=^P.(/)A(/) ■ 

where 

and 7(/) is the overlap reduction function, arising from the different locations and orien- 
tations of the two interferometers. As shown in Figure 2, the identical antenna patterns of 
the collocated Hanford interferometers imply 7(/) = 1. For the Hanford-Livingston pair the 
overlap reduction is significant above 50 Hz. In Equations 5 and 6, Scwif) is the strain 
power spectrum of the stochastic GW background to be searched. Assuming a power-law 
template GW spectrum with index a (see Equation 2), the normalization constant J\f in 
Equation 5 is chosen such that < Y >= Qq,T. 

In order to deal with data non-stationarity, and for purposes of computational feasibility, 
the data for an interferometer pair are divided into many intervals of equal duration, and Yi 
and are calculated for each interval /. The data in each interval are decimated from 16384 
Hz to 1024 Hz and high-pass filtered with a 40 Hz cut-off. They are also Hann- windowed 
to avoid spectral leakage from strong lines present in the data. Since Hann-windowing 
effectively reduces the interval length by 50%, the data intervals are overlapped by 50% to 
recover the original signal-to-noise ratio. The effects of windowing are taken into account as 
discussed in (Abbott et al. 2004b). 

The PSDs for each interval (needed for the calculation of Qi{f) and of ayj) are calculated 
using the two neighboring intervals. This approach avoids a bias that would otherwise exist 
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Fig. 2. — Overlap reduction function for the Hanford-Hanford pair (black solid) and for the 
Hanford- Livingston pair (gray dashed). 
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due to a non-zero covariance between the cross-power Y{f) and the power spectra estimated 
from the same data (Bendat & Piersol 2000). It also allows for a stationarity cut, which we 
will describe in more detail below. 

We consider two interval durations and frequency resolutions: 

• 60-sec duration with 1/4 Hz resolution: the PSDs are calculated by averaging 58 50% 
overlapping periodograms (based on the two neighboring 60-sec intervals) in Welch's 
modified periodogram method. 

• 192-sec duration with 1/32 Hz resolution: the PSDs are calculated by averaging 22 50% 
overlapping periodograms (based on the two neighboring 192-sec intervals) in Welch's 
modified periodogram method. 

As we will discuss below, the 60-sec intervals allow better sensitivity to noise transients 
and are better suited for data-stationarity cuts, while the 192-sec intervals allow higher 
frequency resolution of the power and cross-power spectra and are better suited for removing 
sharp lines from the analysis. 

The data for a given interval / are Fourier transformed and rebinncd to the frequency 
resolution of the optimal filter to complete the calculation of Yj (Eq. 3). Both the PSDs and 
the Fourier transforms of the data arc calibrated using interferometer response functions, 
determined for every minute of data using a measurement of the interferometer response 
to a sinusoidal calibration force. To maximize the signal-to-noise ratio, the intervals are 
combined by performing a weighted average (with weights l/cTy^), properly accounting for 
the 50% overlapping as discussed in (Lazzarini & Romano 2004) . 



The results of this paper are based on the Hanford- Livingston interferometer pairs, 
for which the broadband instrumental correlations are minimized. Nevertheless, it is still 
necessary to investigate if there are any remaining periodic instrumental correlations. We 
do this by calculating the coherence over the whole S4 run. The coherence is defined as: 



2.2. 



Identification of Correlated Instrumental Lines 



r(/) 



Plimif) ■ 



(7) 



The numerator is the square of the cross-spectral density (CSD) between the two interfer- 
ometers, and the denominator contains the two power spectral densities (PSDs). We average 
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Fig. 3. — Coherence calculated for the HlLl pair (top) and for the H2L1 pair (bottom) over 
all of S4 data for 1 mHz rcsohition and 100 mHz resolution. The horizontal dashed lines 
indicate 1/Navg - the expected level of coherence after averaging over N^vg time-periods with 
uncorrelated spectra. The line at 376 Hz is one of the simulated pulsar lines. 
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the CSD and the PSDs over the whole run at two different resolutions: 1 mHz and 100 mHz. 
Figure 3 shows the results of this calculation for the HlLl pair and for the H2L1 pair. 

At 1 mHz resolution, a forest of sharp 1 Hz harmonic lines can be observed. These 
lines were likely caused by the sharp ramp of a one-pulse-per-second signal, injected into the 
data acquisition system to synchronize it with the Global Positioning System (GPS) time 
reference. After the S4 run ended, the sharp ramp signal was replaced by smooth sinusoidal 
signals, with the goal of significantly reducing the 1 Hz harmonic lines in future LIGO data 
runs. In addition to the 1 Hz lines, the 1 mHz coherence plots in Figure 3 also include 
some of the simulated pulsar lines, which were injected into the differential-arm servo of 
the interferometers by physically moving the mirrors. Both the 1 Hz harmonics and the 
simulated pulsar lines can be removed in the final analysis, and we will discuss this further 
in Section 3. Figure 4 shows that the histogram of the coherence at 1 mHz resolution follows 
the expected exponential distribution, if one ignores the 1 Hz harmonics and the simulated 
pulsar lines. 



2.3. Data Quality Cuts 

In our analysis, we include time periods during which both interferometers are in low- 
noise, science mode. We exclude: 

• Time-periods when digitizer signals saturate. 

• 30-sec intervals prior to each lock loss. These intervals are known to be particularly 
noisy. 

We then proceed to calculate Yj and ayj for each interval /, and define three data-quality 
cuts. First, we reject intervals known to contain large glitches in one interferometer. These 
intervals were identified by searching for discontinuities in the PSD trends over the whole 
S4 run. Second, we reject intervals for which ctyj is anomalously large. In particular, for the 
192-sec analysis, we require dy^ < 1 sec for the HlLl pair, and (Tyj < 2 sec for the H2L1 pair 
(recall that Y is normalized such that < Y >= Vt^T ^ with T = 192 sec in this case). The 
glitch cut and the large-sigma cut largely overlap, and arc designed to remove particularly 
noisy time-periods from the analysis. Note, also, that due to the weighting with 1/cry^, the 
contribution of these intervals to the final result would be suppressed, but we reject them from 
the analysis nevertheless. Third, we reject the intervals for which Act = \(Ty, — '7yj\/'7yi > C- 
Here, (7y^ is calculated using the two intervals neighboring interval /, and a'y^ is calculated 
using the interval / itself. The optimization of threshold ( is discussed below. The goal of 
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Fig. 4. — Histogram of the coherence for HlLl (top) and H2L1 (bottom) at 1 mHz resohition 
follows the expected exponential distribution, with exponent coefficient Navg (the number of 
time-periods over which the average is made). 
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this cut is to capture noise-transients in the data, and reject them from the analysis. Figure 
5 shows the impact of these cuts for the HlLl pair, analyzed with 192-sec segments, 1/32 
Hz resolution, and with ( = 0.3. This Figure also shows daily variation in the sensitivity 
to stochastic GW background, arising from the variation in the strain sensitivity of the 
interferometers, which is typically worse during the week-days than during the weekends or 
nights. 

Figure 6 shows the distribution of the residuals for the same analysis. For a given 
interval /, the residual is defined as 

Note that the data quality cuts remove outliers from the residual distribution, hence making 
the data more stationary. After the cuts, the Kolmogorov-Smirnov test indicates that the 
residual distribution is consistent with a Gaussian, for both HlLl and H2L1 analyses with 
192-sec intervals, 1/32 Hz resolution, and ( — 0.3. 



3. Results 

3.1. New Upper Limit 

We performed a "blind" analysis for the HlLl and the H2L1 pairs with 60-sec intervals, 
1/4 Hz resolution, and ( = 0.2. To avoid biasing the results, all data-quality cuts were 
defined based on studies done with a 0.1 sec time-shift between the two interferometers in 
a pair. Such a time-shift removes any GW correlations, without significantly affecting the 
instrumental noise performance. After the data quality cuts were finalized, we made one 
last pass through the data, with zero time-shift, and obtained the final results of the blind 
analysis. 

The results from the blind analysis for the frequency-independent template spectrum 
{a — 0) are listed in the first row of Table 1 for HlLl and in the first row of Table 2 for 
H2L1. These results show no evidence of a stochastic GW background. After completing 
the blind analysis, we discovered that the instrumental 1 Hz harmonic lines, discussed in 
Section 2.2, are correlated between the two sites. We felt compelled on scientific grounds not 
to ignore these correlations, even though they had been discovered after our initial, blind, 
analysis was complete. 

In order to remove from our results any possible infiuence of the correlated lines, we 
repeated our analysis with refined frequency resolution of 1/32 Hz. We increased the interval 
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Fig. 5. — Trend of cry, over the whole S4 run for the 192-sec intervals of HlLl pair. The 
dashed line denotes the large a cut: segments lying above this line are removed from analysis. 
Note the daily variation in the sensitivity of this pair. 




Fig. 6. — Distribution of residuals for the HlLl pair with 192-sec segments: all data are 
shown in gray, data that passes data quality cuts are shown in black, and the Gaussian fit 
to the black histogram is shown as a dashed curve. 
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length from 60-sec to 192-sec, which imphes that the PSDs are estimated by averaging 22 
50% overlapping periodograms. These changes allowed us to exclude the 1 Hz harmonics 
from the analysis while losing only ~ 3% of the bandwidth. The drawback, however, was 
that the 192-sec analysis was less able to identify and exclude the noise-transients than the 
60-sec analysis, as shown in Figure 7. As a result, the Aa cut was retuned for the 192-sec 
analysis after the 60-sec analysis with zero time-shift was completed. 

We examined several methods of defining the data quality cuts for the 192-sec analysis. 
First, we calculated the Aa cut in the 60-sec analysis with ( = 0.2 (along with large sigma 
and glitch cuts), and then declared 192-scc intervals "bad" if they overlapped with the "bad" 
60-sec intervals. This approach has the advantage that the 60-sec analysis is more sensitive 
to transients, and the disadvantage that a significant fraction of the data is rejected. The 
results of this approach are summarized in the second row of Tables 1 and 2. Second, we 
defined all cuts directly on 192-sec intervals. We varied the value of C (0.2, 0.3, and 0.4) and 
selected ( = 0.3 as the optimal parameter choice. This choice represents the best compromise 
between data loss and data quality: it rejects the least amount of data, while still excluding 
all residual outliers and preserving the Gaussianity of the data. The last three rows of Tables 
1 and 2 summarize the results for the three values of 

The 192-sec analysis was not performed blindly. However, it agrees with our blind 
60-sec analysis very well, because the contribution of the correlated 1 Hz harmonic lines is 
only about 1.5a. It is also more conservative than the blind analysis since the value of the 
theoretical error is larger due to a smaller amount of data available in the form of acceptable 
192-sec intervals (as compared to the 60-sec intervals). It also properly handles the known 
instrumental correlations at 1 Hz harmonics. Hence, in the remaining part of the paper we 
will focus on the 192-sec analysis with ( — 0.3. 

Figures 8 and 9 show the results of the 192-sec analysis with ( = 0.3 for HlLl and 
H2L1 respectively. The top plots in these two Figures show the cummulative estimates for 
the frequency-independent GW template {a = 0) as a function of time. They indicate that 
there is no particular time during the run that dominates the result. Moreover, the ±1.65(7 
bounds converge as ~ '^-/VT, as expected. The middle plots of the two Figures show the 
cross-correlation spectra (i.e. the integrand of Equation 3). They indicate that there is no 
particular frequency that dominates the result. 

We perform a weighted average of the HlLl and H2L1 results at each frequency bin (with 
inverse variances as weights). The resulting cross-correlation spectrum and the theoretical 
error are shown in Figure 10. The frequency range 51-150 Hz contributes 99% of the full 
sensitivity, as determined by the inverse combined variance. Integrating over this frequency 
range, we obtain the final estimate for the frequency-independent spectrum: Qq — (—0.8 ± 
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Fig. 7. — Trend of uy for the 60-sec analysis (gray) and 192-sec analysis (black) over a short 
period of time. The two bands were scaled to overlap. Note that the gray band is wider, 
indicating that the 60-sec analysis is more sensitive to noise variations. 
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Fig. 8. — HlLl, 192-sec analysis with ( = 0.3. Top: Cummulative estimate of f2o is shown 
as a function of time. Middle: cross-correlation spectrum Y{f). Bottom: theoretical uncer- 
tainty CTyif) as a function of frequency. 
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Fig. 9. — H2L1, 192-sec analysis with ( = 0.3. Top: Cummulative estimate of f2o is shown 
as a function of time. Middle: cross-correlation spectrum Y{f). Bottom: theoretical uncer- 
tainty CTyif) as a function of frequency. 
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10. — Combined HlLl + H2L1 result, 192-sec analysis with ( ~ 0.3. Top: Combined 
-correlation spectrum. Bottom: Theoretical uncertainty cy(/). 
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50-500 Hz 


51-150 Hz 




60-sec (blind) 


0.2 


7.1% 


6.4 ±4.3 




13.0 


192-sec 


conversion 


21.0% 


3.4 ±4.9 


-2.4 ±5.0 


7.0 


192-sec 


0.2 


10.9% 


6.1 ±4.7 


0.1 ±4.8 


8.0 


192-sec 


0.3 


6.5% 


4.8 ±4.6 


-0.7 ±4.7 


7.3 


192-sec 


0.4 


5.0% 


4.8 ±4.6 


-0.7 ±4.7 


7.3 



Table 1: HlLl results, for hioo = 0.72 and for a frequency-independent template spectrum 
(a = 0). The first two columns define the analysis (interval duration and Act cut). For the 
192-sec analysis with ( denoted as "conversion" (second row), the Act cut was defined using 
the 60-sec analysis listed in the first row: the 192-sec segments are declared "bad" if they 
overlap with a "bad" 60-sec segment. The third column shows the fraction of the data lost 
to data-quality cuts. The fourth column lists the estimates of Qq for the 50-500 Hz range 
before notching the 1 Hz harmonics and the simulated pulsar lines. The fifth column lists 
the estimates of Qq for the 51 — 150 Hz range after notching the 1 Hz harmonics and the 
simulated pulsar lines - this frequency range was determined to be optimal for the combined 
H1L1+H2L1 analysis (see text). The sixth column gives the 90% upper limit based on the 
result in the fifth column. For all analyses presented here, the distribution of the residual 
outliers was consistent with a Gaussian distribution (c.f. Figure 6; the Kolmogorov-Smirnov 
test statistic for comparing the two distributions was greater than 85% in all cases). 
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Table 2: H2L1 results, for /iioo = 0.72 and for a frequency-independent template spectrum 
{a — 0). The columns are as in Table 1. For all analyses presented here, the distribution of 
the residual outliers was consistent with a Gaussian distribution (the Kolmogorov-Smirnov 
test statistic for comparing the two distributions was greater than 52% in all cases). 
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The dominant systematic uncertainty of this result comes from the amphtude cahbra- 

tion uncertainty of the interferometers. This uncertainty is estimated to be 5% (LI) and 8% 
(HI and H2), and it is frequency- independent. The uncertainty in the phase of the interfer- 
ometer strain response is neghgible compared to the magnitude and statistical uncertainties. 
Similarly, the effect of timing errors, measured to be < 4 fis, is negligible. Using hardware 
injections, we estimate that the effect of the timing errors on our point estimate is < 0.2%. 

We then construct the Bayesian posterior distribution for using the above estimate, 
following (Loredo 1990). We assume a Gaussian distribution for the amplitude calibration 
uncertainty (with mean 1 and standard deviation V0.05^ + 0.08^ = 0.093), and wc marginal- 
ize over it. We assume the prior distribution for Qq to be the posterior distribution obtained 
in our previous analysis of the S3 data (Abbott et al. 2005). The 90% upper hmit is the 
value of Jlo for which 90% of the posterior distribution lies between and the upper limit. 
This procedure yields the Bayesian 90% UL on Qq of 6.5 x 10~^. This is an improvement 
by a factor 13 x over the previous result in the same frequency band, established based 
on the science run S3 (Abbott et al. 2005). To investigate robustness of our result under 
different priors, we repeated the calculation using flat priors for the amplitude calibration 
uncertainties, and using a flat prior for Qq between and 8.4 x 10^"^ (previous 90% UL by 
LIGO, (Abbott et al. 2005)). We found that these different choices of priors have less than 
3% effect on the 90% upper hmit. 

Once the estimate is made for the frequency-independent spectrum, one can perform 
appropriate frequency- dependent scalings of Y{f) and (XY{f), recalculate the posterior distri- 
butions, and re-marginalize to obtain upper limits for other templates, such as the power-law 
templates with different spectral indices a (see Equation 2). Figure 11 shows the 90% UL as 
a function of the spectral index a obtained for this analysis. Similar results for the S3 run of 
LIGO, as well as the expected sensitivities of HlLl and H1H2 pairs assuming 1 year of expo- 
sure and design interferometer sensitivities, are also shown. The frequency range of interest 
is deflned to include 99% of the full sensitivity, as determined by the inverse variance. For 
the S4 result, the frequency range varies between (50 — 107) Hz for a = —3 and (73 — 288) 
Hz for a = +3, as shown in the bottom plot of Figure 11. Note that the expected sensitivity 
of the collocated Hanford interferometer pair (H1H2) is signiflcantly better than that of the 
HlLl pair due to the better antenna pattern overlap (see Figure 2). However, this pair is 
also more susceptible to instumentally correlated noise. New analysis methods are being 
pursued to estimate and suppress these instrumental correlations (Fotopoulos 2006). 
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Fig. 11.— Top: 90% UL on as a function of a for S3 HlLl and S4 H1LH-H2L1 com- 
bined, and expected final sensitivities of LIGO HlLl and H1H2 pairs, assuming LIGO design 
sensitivity and one year of exposure. Bottom: Frequency band containing 99% of the full 
sensitivity (as determined by the inverse variance) is plotted as a function of a for the S4 
result. 
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3.2. Signal Injections 

We exercise the analysis procedure described above using simulated stochastic signals 
injected into the data both in software and in hardware - see (Bose et al. 2003). In particular, 
we verify that the recovery of the injected signals is not affected by the data-quality cuts we 
impose. The hardware injections are performed by physically moving the mirrors to simulate 
a stochastic GW signal. Three hardware injections were performed during the S4 run, all 
using a frequency-independent GW spectrum {a — 0). Table 3 summarizes the recovery of all 
of the hardware injections. Figure 12 shows the cross-correlation spectrum for injection 3. It 
also shows the inverse Fourier transform of the spectrum, which is equivalent to the estimate 
of flo for different values of time-lag between two interferometers (for short time-lags). 



Injection 


HlLl Expected 


HlLl Recovered 


H2L1 Expected 


H2L1 Recovered 




(xlO-2) 


(xlO-2) 


(xlO-2) 


(xlO-2) 


1 


9.1 


7.9 ±0.2 ±0.8 


7.9 


6.9 ±0.4 ±0.6 


2 


2.5 


2.5 ±0.4 ±0.2 


2.3 


1.5 ±0.6 ±0.2 


3 


1.1 


0.95 ±0.04 ±0.10 







Table 3: Summary of the hardware injection amplitudes during S4. The second and fourth 
columns indicate the expected injection amplitudes for HlLl and H2L1 respectively, based 
on the signal injected into the differential-arm servo. The third and the fifth columns list 
the recovered values using the HlLl and H2L1 pairs respectively. The recovered values are 
listed with statistical errors (as defined in Equation 4), and with systematic errors (estimated 
using 5% calibration uncertainty in LI and 8% calibration uncertainty in HI and H2, added 
in quadrature). For injection 3, the data of H2 interferometer were compromised due to a 
failure of the interferometer's laser. 

We performed a sequence of software injections, where the injected signal is simply added 
to the interferometer data in the analysis. We performed 10 trials for 4 injection amplitudes 
using about 1/3 of the S4 HlLl data. Figure 13 shows that the signal is successfully recovered 
down to 1 X 10~^ in Qq. Moreover, the theoretical error bars agree well with the standard 
deviation over the 10 trials. All injections were performed assuming a frequency-independent 
GW spectrum {a — 0). 

4. Implications 

In this Section, we investigate the implications of the new upper limit for some of the 
models of the stochastic GW background. We also discuss the complementarity of our result 
with various other experimental constraints on the stochastic GW background. 
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Fig. 12. — HlLl injection 3 (with intended injection amplitude Qq = 1.1 x 10~^). Top: 
Cross-correlation spectrum Y{f). Bottom: Inverse Fourier transform of the cross-correlation 
spectrum indicates the clear signal at zero lag. The inset plot is a zoom-in around zero lag. 
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Fig. 13. — Injections using HlLl data: 10 trials were performed for software injections with 
amplitudes = 1 x 10^^, 2 x 10^"^, 6 x lO^"', and 2 x 10^^. The left(gray) error bars denote 
the theoretical errors, while the right (black) error bars denote the standard errors over the 
10 trials. The remaining points correspond to the three hardware injections listed in Table 
3; their error bars correspond to statistical and systematic errors added in quadrature, as 
shown in Table 3. 
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4.1. Complementarity with Other Measurements and Observations 

Figure 14 compares different experiments and some of tlie tlieoretical models. For 
wavelengtlis larger than the horizon size at the surface of last scattering (redshifted to today, 
this corresponds to frequencies below ~ 10^^^ Hz), the COBE observations of the CMB place 
an upper limit on the stochastic gravitational wave background of flcwif) < 1-3 x 10~^^ 
(Allen & Koranda 1994). In standard inflationary models (Turner 1997), the GW spectrum 
is expected to be (almost) flat at frequencies above ~ 10~^^ Hz. 

The fluctuations in the arrival times of millisecond pulsar signals can be used to place 
a bound at ~ 10~^ Hz (Jenet et al. 2006): ^Gw{f) < 3.9 x 10~^ (assuming frequency 
independent GW spectrum). Similarly, Doppler tracking of the Cassini spacecraft can be 
used to arrive at yet another bound, in the 10~^ — 10~^ Hz band (Armstrong et al. 2003): 



If the energy density carried by the gravitational waves at the time of Big-Bang Nucle- 
osynthesis (BBN) were large, the amounts of the hght nuclei produced in the process could 
be altered. Hence, the BBN model and observations can be used to constrain the total 
energy carried by gravitational waves at the time of nucleosynthesis (Kolb & Turner 1990; 
Maggiore 2000; Allen 1996): 



where N,, is the effective number of relativistic species at the time of BBN. Measurements of 
the light-element abundances, combined with the WMAP data, give the following 95% upper 
bound: A'^^ — 3 < 1.4 (Cyburt et al. 2005). This limit translates into J ficwl/) d{\nf ) < 
1.5 X 10"^. This bound applies down to ~ 10~^° Hz, corresponding to the horizon size at 
the time of BBN. 

Gravitational waves are also expected to leave a possible imprint on the CMB and matter 
spectra, similar to that of massless neutrinos. (Smith et al. 2006a) used recent measurements 
of the CMB anisotropy spectrum, galaxy power spectrum, and of the Lyman— a forest, to 
constrain the energy density carried by gravitational waves to J ^owif) d{ln.f) < 1.3 x 10~^ 
for homogeneous initial conditions. This bound is competitive with the BBN bound and 
it extends down to ~ 10~^^ Hz, corresponding to the horizon size at the time of CMB 
decoupling. It is also expected to improve as new experiments come online (such as Planck 
or CMBPol). 

The LIGO results apply to the frequency region around 100 Hz. The result discussed in 
this paper is an improvement by a factor 13 x over the previous LIGO result in the 100 Hz 
region, for a frequency-independent spectrum of GW background. A one-year run at design 



^^Gw(/) < 0.027. 
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Fig. 14. — Landscape plot (see text for details). The curves corresponding to inflationary, 
cosmic-string, and pre-big-bang models are examples; significant variations of the predicted 

spectra are possible as the model parameters are varied. The bounds labeled "BBN" and 
"CMB and Matter Spectra" apply to the integral of the GW spectrum over the frequency 
range spanned by the corresponding lines. 
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sensitivity of LIGO (the S5 run, which began in November of 2005) is expected to improve 
the sensitivity by another factor 10 x - 100 x, while Advanced LIGO is expected to achieve 
sensitivities better by yet another factor lOOx - lOOOx, eventually reaching 10"^^ — 10^^ for 
Qq. The uncertainty in the final reach of LIGO and of Advanced LIGO comes from the 
potential instrumental correlations that could be present between the collocated Hanford 
interferometers and from the uncertainty in the final optical configurations of interferometers 
in Advanced LIGO. The result discussed in this paper is still weaker than the indirect BBN 
bound, but the future runs by LIGO and Advanced LIGO are expected to surpass this 
bound. 

The standard inflationary models are most stringently constrained by the CMB bound 
at lowest frequencies. Although they are most likely out of range of LIGO and of Advanced 
LIGO, they may be accessible to future GW interferometers (Smith et al. 2006b). However, 
there are models of stochastic GW background that LIGO is beginning to explore. We 
illustrate this with examples of cosmic string and pre-big-bang models. 

4.2. Implications for Cosmic Strings Models 

Cosmic strings can be formed as linear defects during symmetry breaking phase transi- 
tions in the early Universe as well as in string theory inspired inflation scenarios. In the latter 
case they have been dubbed cosmic superstrings. CMB data is not consistent with cosmic 
strings as the predominant source of density fluctuations in the Universe. Their existence, 
however, is not ruled out below the GUT scale and cosmic strings may still lead to a myriad 
of detectable astrophysical signatures such as gravitational radiation, gamma-ray bursts and 
ultra-high energy cosmic rays. For a review, see (Vilenkin & Shellard 2000). 

(Damour & Vilenkin 2005) investigated the stochastic background of gravitational waves 
produced by cusps on cosmic strings, integrated over all redshifts and all directions on the 
sky. They flnd that three parameters deflne the gravitational-wave spectrum due to cosmic 
strings: 

• String tension /i: This parameter is usually expressed as a dimensionless quantity G/i 
(assuming speed of light c—1), where G is Newton's constant. String-theory inspired 
inflation scenarios prefer the range 10~^^ < G/j, < 10~^. 

• Reconection probability p: While ordinary, fleld-theoretic strings reconnect when- 
ever they intersect (p = 1), the reconnection probability for superstrings is typi- 
cally smaller than 1. In particular, the theoretically favored range is 10~^ < p < 1 
(Jackson et al. 2005). 
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• e: This parameter describes the typical size of the closed loops produced in the string 
network. The value of this parameter is uncertain, and it can span several orders of 
magnitude. We will consider the range 10~^^ < e < 10~^, which is both theoretically 
viable and most interesting from the point of view of LIGO. 

These parameters determine both the amplitude and the shape of the gravitational- 
wave spectrum. In particular, parameters e and G/x determine the lowest frequency (at a 
given redshift) at which a string loop could emit gravitational radiation. Since there is a 
low-frequency cut-off to the predicted gravitational-wave spectrum, it is possible to have 
a cosmic string model that would avoid the low-frequency bounds due to CMB or pulsar 
timing measurements, but still be within reach of LIGO (c.f. Figure 14). 

Figure 15 shows the region of the parameter space (for p = 10~^) excluded by the result 
discussed in this paper and by the S3 result (Abbott et al. 2005), as well as the expected 
reach of LIGO and Advanced LIGO in future runs. As shown in this Figure, LIGO is most 
sensitive to the regions of large Gfi and small e. Moreover, LIGO results are complementary 
to the pulsar timing limit, which is most sensitive to models with large Gn and large e, and to 
the indirect BBN limit. In particular, the population of models with p = 10^'^, e ^ 5 x 10^^^, 
and 5 x 10~^ < On ^ 7 x 10~^ are excluded by the result discussed in this paper, but are not 
accessible to other current experimental bounds. To produce this Figure, we used Equations 
4.1-4.7 from (Damour & Vilenkin 2005). 

We note that as p increases toward 1, the bounds/reach of all experiments in the 
Gfi — e plane weaken because the amplitude of the gravitational-wave spectrum scales with 
1/p. We also note that there are significant uncertainties in the calculation discussed by 
(Damour & Vilenkin 2005), mostly due to incomplete understanding of the string network 
behavior. For example, the size of loops is usually assumed to be given by the scale 
of gravitational back reaction. Recent numerical simulations of cosmic string networks 
(Martins & Shellard 2006), (Vanchurin et al. 2005), and (Ringeval et al. 2005) instead sug- 
gest that loop sizes are related to the large scale dynamics of the network. If this is the case 
then the expected stochastic background at the frequencies of pulsar timing experiments 
could be subtantially larger (Hogan 2006). Additionally, the treatment of (Damour & Vilenkin 2005) 
did not account for the effects of late-time acceleration. Thus a more detailed analysis of 
the stochastic background produced by cosmic strings is necessary. 
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p= 10" 




Fig. 15. — The e - Gji plane for the cosmic string models with p — 10~^. The shaded 
regions are excluded by the LIGO S3 limit (darker) and by the LIGO S4 hmit presented 
here (lighter). The hatched regions are accessible to future LIGO runs: expected LIGO 
sensitivity for the HlLl pair, assuming design interferometer strain sensitivity, and 1 year of 
exposure ('—'); expected LIGO sensitivity for the H1H2 pair, assuming design interferometer 
strain sensitivity, and 1 year of exposure ('\'); expected Advanced LIGO sensitivity for the 
HI 112 pair, assuming interferometer strain sensitivity tuned for the binary neutron star 
inspiral search and 1 year of exposure ('/')• The dash-dotted black curve is the exclusion 
curve based on the pulsar limit (Jenet et al. 2006) (the excluded region is above the curve). 
The solid black curve is the exclusion curve based on the indirect big-bang-nucleosynthesis 
bound (the excluded region is to the right of the curve) . 
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4.3. Implications for Pre-Big-Bang Models 

Pre-big-bang models are cosmology models motivated by string theory (Gasperini & Veneziano 1993), 
(Gasperini & Veneziano 2003). In these models, the Universe evolves through several phases: 
"dilaton" phase in which the Universe is large and shrinking; "stringy" phase in which the 
curvature of the Universe is high; and standard radiation and matter dominated phases. The 
GW spectrum is generated by amplification of vacuum fiuctuations as the Universe transi- 
tions from one phase to another. The shape and amplitude of the spectrum are determined 
by the states of the Universe in the different phases. Although the "stringy phase" of the 
model and the transition to the radiation phase are not well understood, some models have 
been proposed in the literature that may partially describe it. In the formalism developed 
by (Buonanno et al. 1997), the GW spectrum produced by the model can be described as: 

• ficwif) ~ P for / < fs, where fg is essentially unconstrained. 

• f^Gw(/) ~ y-^^^AJ £qj. fs < f < fi, where // < 1.5 defines the evolution of the Universe in 
the "stringy" phase. The cutoff frequency /i is determined by string-related parameters 
and its most natural value is expected to be 4.3 x 10^° Hz. 

Following the analysis in (Mandic & Buonanno 2006), we scan the parameter space 
(/s, /X, /i) . For each set of values of these parameters, we calculate the GW spectrum following 
(Buonanno et al. 1997) and check whether the model is accessible to the current (and past) 
LIGO results. We also project the sensitivity of initial and Advanced LIGO to these models. 
Figure 16 shows the /i — plane for fg = 30 Hz. Note that the LIGO S3 and S4 results are 
beginning to explore the parameter space of these models, although the indirect BBN bound 
is still a stronger constraint. Future runs of LIGO and of Advanced LIGO are expected 
to explore significantly larger parts of the parameter space, eventually surpassing the BBN 
bound (in some parts of the parameter space) and even reaching the most natural value of 
/i = 4.3 X 10^° Hz. 

5. Conclusions 

LIGO data acquired during the science run S4 yield a new Bayesian 90% upper limit 
on the amplitude of the stochastic GW background: Qq < 6.5 x 10~^ for the frequency- 
independent GW spectrum {a — 0) in the frequency band 51-150 Hz. Similar limits are 
obtained for other values of a, as shown in Figure 11. This result is an improvement by a 
factor 13 X over the previous upper limit in the same frequency range, obtained by LIGO in 
the science run S3. 



10 

1,2 1,25 1,3 1,35 1,4 145 1,5 

Fig. 16. — The /i — fi plane for the pre-big-bang models with fs = 30 Hz. The shaded 
regions are excluded by the LIGO S3 limit (darker) and by the LIGO S4 limit presented 
here (lighter). The hatched regions are accessible to future LIGO runs: expected LIGO 
sensitivity for the HlLl pair, assuming design interferometer strain sensitivity, and 1 year of 
exposure ('—'); expected LIGO sensitivity for the H1H2 pair, assuming design interferometer 
strain sensitivity, and 1 year of exposure ('\'); expected Advanced LIGO sensitivity for the 
H1H2 pair, assuming interferometer strain sensitivity tuned for the binary neutron star 
inspiral search and 1 year of exposure ('/'). The solid black curve is the exclusion curve 
based on the BBN limit (the excluded region is above the curve). The horizontal dashed line 
denotes the most natural value of /i = 4.3 x 10^° Hz. 
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This result is obtained using 192-sec long intervals of data with 1/32 Hz frequency 
resolution of the spectra, and it properly excludes the known instrumental correlations at 1 
Hz harmonics. It is fully consistent with the blind result that uses 60-sec long intervals of 
data with 1/4 Hz frequency resolution, which is slightly contaminated by the instrumental 
1 Hz harmonics. It is also more conservative than the blind result, as the theoretical error is 
larger due to the smaller amount of data available in the form of acceptable 192-sec intervals 
(as compared to the 60-sec intervals). 

This result is complementary to the constraints on the gravitational-wave spectrum, 
based on the measurements of the CMB spectrum and pulsar timing. It is still weaker 
than the indirect BBN bound in the relevant frequency range. The ongoing 1-year long 
run of LIGO at the design sensitivity, and the future runs of Advanced LIGO, are expected 
to surpass the BBN bound. Furthermore, this result is already exploring the parameter 
space of some models of the stochastic GW background, such as cosmic strings models and 
pre-big-bang models. 
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